Selective coherent x-ray diffractive imaging of displacement fields in 
(Ga,Mn)As/GaAs periodical wires 
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Coherent x-ray diffractive imaging is extended to high resolution strain analysis in crystalline nanos- 
tructured devices. The application potential is demonstrated by determining the strain distribution 
in (Ga,Mn)As/GaAs nanowires. By separating diffraction signals in reciprocal spaces, individual 
parts of the device could be reconstructed independently by our inversion procedure. We demon- 
strate the method to be effective for material specific reconstruction of strain distribution in highly 
integrated devices. 

PACS numbers: 61.05.cp, 62.20.-x, 42.30.Rx, 75.50.Pp 



Nanoscopic elastic strain analysis is crucial for under- 
standing physical properties and fabrication processes 
of crystalline nanostructures and devices. In epitaxial 
nanostructures, strain distribution is strongly correlated 
to their shape and composition. It is also an important 
driving force for self-assembling processes 

In conventional x-ray diffraction experiments informa- 
tion about shape, composition and strain is obtained 
indirectly by quantitative comparison of measured and 
simulated reciprocal-space distributions of scattered in- 
tensity 0. Simulations are based on a model with pre- 
liminary assumptions about detailed shape, composition 
and strain profile, which will be confirmed by the exper- 
imental results as long as the model structure fits to the 
reality. 

A direct determination of the density of the scatter- 
ing centers from the reciprocal-space distribution of the 
scattered intensity is however not possible, unless a spe- 
cial phase- retrieval method is used. This phase-retrieval 
procedure known as coherent x-ray diffractive imaging 
(CXDI) has a potential to become a powerful technique 
for the structural characterization of small objects and 
micro-sized semiconductor devices 0, 0, 0, 0, 0,M S Efll • 
The method principle is based on the iterative loop of 
direct and inverse Fourier transformation (FT) (towards 
the experimental intensity distribution and back to the 
sample space) and it may refine the genuine object den- 
sity distribution g(r) even by starting from a model 



which is far from reality [111. Il2|. 



In case of crystalline nano-objects the distribution of 
the scattered amplitude in reciprocal space around a re- 
ciprocal lattice point (RLP) is given by 
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Here, h and Q denote the reciprocal-lattice vector and 
the position vector in reciprocal space (scattering vector), 
respectively, Qh(r) is the corresponding Fourier compo- 
nent of electron density and u(r) is the displacement field 
with respect to the non-strained lattice. 

Once the phase of the scattered radiation is known, 
direct space amplitude Qh(r) and phase h.u(r) can be 
obtained by inverse FT, within the validity of the kine- 
matical scattering theory (first Born approximation). 

An important step in the phase-retrieval loop is the 
application of the direct (real)-space constraints (see [Uj, 
among others). Usually, the direct-space constraint has 
the form of an a-priori known support function determin- 
ing the region in direct space, where the complex electron 
density g(r) = Qh(?) exp(— \h.u(r)) may differ from zero. 
The first demonstration of the possibility of strain imag- 
ing by phase retrieval using support constraint on par- 
ticular example of weakly strained Pb nanocrystal was 
recently illustrated in [7]. Later it has been shown, how- 
ever, that the knowledge of the support constraint alone 
is not generally sufficient for the phase retrieval of diffrac- 
tion pattern from highly strained crystals [13], but the 
problem may be overcome by introducing additional con- 
straints. Such a suitable constraint in direct space is the 
use of a-priori knowledge to limit the allowed values of 
strain in the iterative process 0, [l3[ . 

The diffraction amplitudes of compositional layered 
nanoobjects can be understood as an interference sum 
of individual wave field contributions scattered from the 
different material regions. We may define firstly a non- 
perturbed state of the nanoobject by assuming for each 
material region idealized homogeneous composition and 
strain, spatially limited by the shapes of the individual 
nanostructures. Secondly we define the perturbation by 
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FIG. 1: Schematic diagram of the iterative loop of our phase 
retrieval algorithm. r and r' in " direct (real)-space con- 
straints" box correspond to the neighbouring points within 
one constituent crystal and e defines the threshold of electron 
density uniformity. i p and i' p are the neighbouring points 
along p-direction (normally lateral or vertical) and Ap corre- 
sponds to the step along this direction (achieved resolution). 



taking now the effects of non-uniform strain relaxation 
and possible compositional fluctuations within previously 
homogeneous layer region of the crystals into account. 
Then, each wave field contribution is the result of con- 
volution of two FTs: One is the individual reciprocal 
shape function which leads to intensity decreases at least 
with q~ 2 with respect to the unperturbed RLP of that 
material. The other one corresponds to the FT of g(r) 
containing the spatial fluctuations in amplitude (compo- 
sitional fluctuation) and phase (displacement field). The 
spread of the later FT in frequency space is dependent 
on the maximum gradient values of displacements in the 
layers and related to the strain-induced broadening of 
diffracted signal AQ around RLP (AQ P - along par- 
ticular p direction). That fact allows to transform the 
"limited frequency constraint" into the direct space con- 
straint of limited component of the displacement gradi- 
ent. Therefore, in addition to the standard support con- 
straint we implement complementary direct-space con- 
straints, namely, the uniformity of electron density (am- 
plitudes) for each constituent crystal component and the 
requirement of continuity of the displacement variations 
(of phase) (see Fig. [1]). 

Now, let us consider the constraint for the vertical fre- 
quencies. In elastically strained-layer nano-objects the 
layer regions of different composition have considerably 
different averaged vertical lattice parameters. The re- 
sulting mean vertical lattice misfit between these regions 
leads to the separation of their diffraction maximums in 
^-direction (see Fig. [3]) and to a shift between their whole 
individual contributions in reciprocal space by Ah. 

Mostly Ah is sufficiently large compared to the maxi- 
mum gradient of displacement field within the individual 
layers. If it is also sufficiently large compared to the main 
contribution of the reciprocal shape functions, recipro- 
cal space regions are distinguishable, where the major 
influence of scattered intensity stems mainly from one 
individual material type. That consequently permits to 
reconstruct g for the individual constituent parts of the 
nanoobject separately. 




FIG. 2: (Color online) Schematic cross-section of the GaM- 
nAs/GaAs periodical wires perpendicular to their direction 
(in the plane of the diffraction experiment). The direct space 
(DS) area attained by the periodical object diffracted signal 
sampling in RS is marked by the dashed line region. 



Experimentally, CXDI can be applied to individual 
nanoobjects like wires and dots, if sufficient high photon 
flux density is available. In order to gain in the radiation 
flux density, microfocusing optics could be used [14] . An- 
other way to overcome the problem is the use of period- 
ical objects. If the scattering coherent volume contains 
a large number of equally spaced identical objects, we 
gain constructive interference concentrating and ampli- 
fying diffracted intensity into so-called grating rods [2[. 
These rods are separated by Aq x = 27r/T, where T is a 
period size. The maximum intensity of the scattered radi- 
ation in rods is given by 7 exp = Ji-^-AT 2 = IiNN c where 
N is the number of the periods irradiated by the primary 
wave, N c is the number of coherently irradiated periods, 
and 1\ is the maximum intensity scattered from a single 
wire. Therefore, if we use an ideally coherent primary 
wave, the scattered intensity is substantially enhanced, 
since it is proportional to 2 . Consequently sampling in 
reciprocal space can be performed in a wider range, where 
measured intensities are without coherent magnification 
already negligibly small. 

Strain analysis in III-V semiconductor nanostruc- 
tures composed of GaAs and (Ga,Mn)As is performed 
to demonstrate the potential of CXDI for invetigation 
of application-relevant nanoobjects. Diluted magnetic 
semiconductors, of which (Ga,Mn)As is a prominent ex- 
ample, have acquired great attention for tunability of 
their magnetic properties and thus for their potential 
in spintronic applications [HI, [l6[ • It has been shown 
that, e.g., their magnetocrystalline anisotropy can be 
substantially affected by elastic strain [17]. One way to 
induce controlled spatial variation of the strain is litho- 
graphic patterning of the GaMnAs layer after the epitax- 
ial growth. Due to the positive lattice mismatch between 
(Ga,Mn)As and GaAs, depending mainly on the concen- 
tration of Mn interstitial atoms [18[, epitaxial (Ga,Mn)As 
layers deposited on GaAs substrates are biaxially com- 
pressed. After creation of a periodic sequence of thin 
wires and groves (surface grating), the film can relax 
freely at the edges of the wires. 

In the investigated structures the period of wires is 
T = 2 fim and the nominal wire width is S x = 1 /im. 
A 200 nm thick (Ga,Mn)As layer with a nominal Mn 
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FIG. 3: (Color online) Measured reciprocal space map (log 
10 of intensities) from (Ga,Mn)As/GaAs periodical wires. Qs 
and Ql confine the diffuse scattering regions from GaAs and 
(Ga,Mn)As respectively. 

concentration of about 7% was grown onto a (001) GaAs 
substrate. The periodic wire structure was prepared by 
electron-beam lithography and reactive ion etching with 
an etching depth of approx. 700 nm (see Fig. [2|). The 
(Ga,Mn)As surface wires were oriented along the [110] di- 
rection. The technological details can be found elsewhere 

The wires do not contain large structure defects like 
dislocations, stacking faults, or precipitates that would 
give rise to considerable x-ray diffuse scattering. X-ray 
diffuse scattering from point-like defects such as Aspa 
antisites, Mn substitutional and interstitial atoms [18[ 
lies far below the detection limit of the experimental set- 
up. Therefore, from the view-point of the experimental 
method used, the wires represent an ideal crystal struc- 
ture with a two-dimensional distribution of the elastic 
displacement vector u(x,z) in the plane xz perpendicu- 
lar to the wire direction. 

Diffraction measurements were performed at ID10B 
beamline at the ESRF. The x-rays from an in- vacuum un- 
dulator source were monochromatized by a diamond(lll) 
crystal to E = 8 keV. The cross-section of the primary 
beam was 0.2x0.2 mm 2 . High resolution reciprocal-space 
images (RSI) were measured by use of a Si(lll) crystal 
analyzer. The result in vicinity of the 004 RLP is shown 
in Fig. [3j The RSI q x q z plane is perpendicular to the 
alignment of the wires. 

Two well separated features in the 004 RSI correspond 
to the major influence areas of the 004 reciprocal lattice 
points of GaAs and (Ga,Mn)As respectively (the corre- 
sponding difference in Bragg angles is about 0.2 degree). 
Grating maxima exhibiting a horizontal spacing of about 
~ 3*10 _4 A _1 can be clearly resolved proving the in-plane 
coherence length of the incoming x-ray beam to be much 



larger than the 2 fim wire period (see the magnified inset 
in Fig. [3]). Since the lower sampling limit of RSM enve- 
lope corresponds to the rod periodicity, the oversampling 
ratio of the measured signal in q x direction is defined as 
&x — • In our case the wire width S x occupies the half 
of the grating period T and consequently the oversam- 
pling rate is about cr x = 2 in q x direction. The internal 
symmetry of the wire in lateral direction decreases the 
unknown information twice. 

The dynamical diffraction peak of the semi-infinite 
GaAs substrate cannot be taken into account since this 
region violates the presumption of the phase retrieval 
method. The eliminated experimental data near the 
GaAs reciprocal lattice point (marked by fl pe ak in Fi g. Op 
allowed to vary freely during reconstruction process [19[ . 

The support in direct space is chosen on the basis of 
growth parameters and is finalized using Parseval theo- 
rem being applied separately to the scattering regions l 
and ^5 in RSI shown in Fig. 03 Performing the phase- 
retrieval algorithm the Fouier components of electron 
densities ph(r) of (Ga,Mn)As and GaAs are assumed to 
be not affected by elastic deformation and shall vary only 
with composition. 

From the diffraction pattern one can conclude on a 
large vertical misfit between the (Ga,Mn)As and GaAs 
regions, allowing us to divide the measured RSI in two 
regions £Il and Vts (see Fig.[3j). Each of them serves fur- 
ther for independent reconstruction of phase and ampli- 
tudes in the corresponding composite parts (Ga,Mn)As 
and GaAs of the wire, respectively. The results are de- 
picted in Fig. [U The reproducibility of the obtained 
solutions was proven performing a series of inversion 
cycles starting from different sets of random phases of 
the scattered radiation. The error metric is defined 

as E l = E?=i(\K calc \ - \/W^?IY^=i I i ieas i where 

\F^ alc \ is the magnitude of the calculated amplitude and 
jmeas - ls ^ e measured intensity of point i in the RSI. 
The error metric is found to be on the level of 10 -3 for 
the (Ga,Mn)As part and becomes even smaller, around 
10 -6 , when inverting the GaAs part. The direct-space 
resolution achieved is about 10 x 10 nm 2 . 

From the results of the independent phase reconstruc- 
tions of the (Ga,Mn)As and GaAs volumes of the wires 
the following observations can be made: 

(i) During the reconstruction of the GaAs part the av- 
erage amplitude in the free-standing wire appears to be 
close to that one in the (Ga,Mn)As (see Fig. [4j). This 
is in good agreement with the fact that the amplitudes 
of structure factors of GaAs and (Ga,Mn)As unit cells 
containing 7% Mn differ by less than 1 percent. 

(ii) The variation of the reconstructed phase fields (i.e., 
the fields of the relative vertical displacement component 
Fig. E](a)) behave continuous at the boundary between 
the (Ga,Mn)As and GaAs volumes. That fact further 
supports the validity of our reconstruction process. 

(iii) Finally, the combination of both independent re- 




FIG. 4: (Color online) Results of reconstructions (a) of Ql 
area of the RSI (see Fig. [3} corresponding to the (Ga,Mn)As 
part, (b) of Qs area corresponding to the GaAs part. For the 
dimensions of direct space of periodic objects see Fig. [2] 




FIG. 5: (Color online) Contour maps of the vertical compo- 
nent u z of the displacement vector in the (Ga,Mn)As/GaAs 
wire relative to the wire state before the relaxation (in A) (a) 
directly reconstructed (calculated from the two phases maps 
shown in Fig. [4j, (b) obtained numerically by finite-element 
simulations. 



constructions into one leads after taking the FT to the 
same RSI as in the Fig. [3] down to the error of 3 * 10 -2 . 
This shows that even overlapping regions in RSI can be 
independently reconstructed, when the difference in in- 
tensity levels between the intensity maxima and the edge 
of the cut map is at least 2 orders of magnitudes. 

We do not interpret the non-homogeneous amplitudes 
in the GaAs substrate part (under the wire). For this 
region strain variation and compositional variation are 
certainly very small. Therefore the missing signal data 
in the RSI close to the strong dynamical peak (excluded 



from the diffraction pattern to be reconstructed) is an 
essential loss of information and leads to expected insta- 
bilities in the reconstructed amplitudes in Fig. H(b). 

Numerical finite-element method (FEM) calulations 
were performed to simulate the strain field inside the 
(Ga,Mn)As/GaAs wires. Simulation used the experi- 
mentally verified value of 0.35% for the lattice mismatch 
between the GaAs substrate and the (Ga,Mn)As layer 
containing 7% Mn. The Poisson ratio and the Young 
modulus of GaAs were used as an approximation for 
(Ga,Mn)As. The results of the simulations are shown 
in Fig. [Sfb). The excellent agreement between simulated 
(FEM) and reconstructed strain field distributions proves 
the high quality of the reconstruction process. 

In summary, we illustrated the ability of CXDI based 
on our phase-retrieval algorithm to study non-uniform 
strain relaxation in strained layer nanoobjects and for 
the first time succeeded to reconstruct the related dis- 
placement fields. Reproducible convergence of the CXDI 
method has been obtained even for strongly inhomoge- 
neously strained objects. The reconstruction procedure 
takes benefit from the material-specific nature of CXDI 
arising from compositionally driven lattice misfit. That 
enables individual analysis of the corresponding compo- 
sitional constituents of the structure. The continuous 
field of displacement can be also recovered by building 
up the strained inhomogeneous system as a whole. The 
quantitative values of displacement field reconstruction 
have been validated using the elasticity theory. Sub ten 
nanometer resolution of the displacement field in direct 
space was achieved. The availability of a robust direct re- 
construction method will allow to gain precise knowledge 
about multicompositional strained nanoobjects. 
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